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Abstract 

We prove the global existence and uniqueness of the classical (weak) solution for the 2D or 
3D compressible Navier-Stokes equations with a density-dependent viscosity coefficient (A = A(p)). 
Initial data and solutions are only small in the energy-norm. We also give a description of the large 
time behavior of the solution. Then, we study the propagation of singularities in solutions. We obtain 
that if there is a vacuum domain at initially, then the vacuum domain will exists for all time, and 
vanishes as time goes to infinity. 

1 Introduction 

In this paper, we consider the following compressible Navier-Stokes equations 

pt + div(pu) = 0, , . 

{pu)t + div(pu (g) u) + VP = ^Am + V((/i + A(p))divu) + pf, ^ ' ' 

for X e and i > 0, = 2 or 3, with the boundary and initial conditions 

u{x,t) — > 0, p{x,t) p> 0, as |x| — > cxD, t > 0, (1.2) 

{p,u)\t=o = {po,uo). (1.3) 

Here p{x,t), u{x,t) and P = P{p) stand for the fluid density, velocity and pressure respectively, / is 
a given external force, the dynamic viscosity coefficient /i is a positive constant, the second viscosity 
coefficient X = \{p) is a function of p. 

In [20] , we proved the global existence of weak solutions for the two-dimensional system, and study the 
propagation of singularities in solutions. In this paper, we want to obtain the global existence, uniqueness 
and the large time behavior of the classical solution to the system (ll.ll) - (|1.3l) in or M^, also obtain 
the global existence of weak solutions and study the propagation of singularities in solutions in R^. 

At first, we obtain the global existence, uniqueness and the large time behavior of the classical 
solution, when the energy of initial data is small, but the oscillation is arbitrarily large. Specifically, we 
fix a positive constant p, assume that (po — P,uo) are small in L^, and pQ — p,uo e with no restrictions 
on their norms, (since we use the classical analysis methods in this paper, we restrict the result of the 
existence of the classical solutions on the framework of Hilbert space H^{]K^) ^ C^(IR^)). Our existence 
result accommodates a wide class of pressures P, including pressures that are not monotone in p. It also 
generalizes and improves upon earlier results of Danchin 4 and Matsumura-Nishida |13| in a significant 
way: (po — p, uq) are only small in L^. 

Now, we give a precise formulation of our result. Concerning the pressure P, viscosity coefficients p, 
and A, we fix < p < p and assume that 

P e C\[0,p]) n C3((0,p]), A e ^^([O, p]) n 0^(0, p]), 

„ , f [0,oo), = 2, , ,^ 
Ae| |q^3^^^^ ^^3^ forallpe[0,p], 

P(0) = 0, P'(p) > 0, 

{p-p)[P{p)-P{p)]>0, pe[0,p)u {p, p] . 
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Let G be the potential energy density, defined by 

Then, for any g G C^([0, p]) with g{p) — g'{p) — 0, there is a constant C such that \g{p)\ < CG{p), 
P&[0,p]. 
Define 

Co = J QpokoP + G(po)) da., (1.6) 
Cf = sup||/(.,i)|li. + / {\\.f{;t)U. + ||/(-,t)|li. +a=^^-^|lV/|li. +a2^-i||/,(.,t)|li.) dt (1.7) 

and 



Mq 



(<^'\\ft\\h + ^'W^fWU + 'y"' ll/tllit^, + cjP' ||V/||i^,) dt + llV^olli. + sup (1.8) 

^ ^ t>0 

where a{t) — min{l,<}, 

r 2 + f, / 4 + 2q, iV = 2, 

^^ = 1 1 + 1, ^ = 3, (1-9) 

and (7 is a constant satisfying 

and,^<-J^,Vpe[0,p-]. (1.10) 

As in [6], we recall the definition of the vorticity matrix w^'*' — dku^ — dju^ , and definition of the 
function 

F= (A + 2/i)divM-P(p) + P(p). (1.11) 

Thus, we have 

= ^.( ^^+"2^^^ ^ + ^^-^^^ 
We also define the convective derivative by — w ~ wt + u ■ \/w, the Holder norm 

^ |i^(a:) - v{y)\ 

< « >A= sup — ^ 

x,yeA F ~ y\ 

and 

a,/3 _ \gix,t)-giy,s)\ 

(x,t)7^(y,s) 

where v.ACR'^^ R^, g : A x [ii,i2] ^ and a,/3 G (0, 1]. 
The following is the main result of this paper. 

Theorem 1.1. Assume that conditions ^1.4^ - (T7^) hold. Then, for a given positive number M (not 
necessarily small) and pi G {p,p), there are positive number e, such that, the Cauchy problem lll.l\) ~ nT3\) 
with the initial data (pq^uq) and external force f satisfying 



< < inf pq < suppo < Pi, 

Co + Cf < e, , . „s 

M, < M, ^^-^^^ 
^ po-p, uoeH^ e C([0,oo);i2), /e C([0,oo);ij2), 
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has a unique global classical solution (p, u) satisfying 

1 



^P, < Pix,t) < p, x€W\ t>0, (1.14) 
{p-p,u) eC^R'^ x[0,T])nC{[0,T],H^)nC\[0,T],H^), (1.15) 

snp i\\ip-p,u)\\H^ + \\{puUt)\\H^)+ [ \\u\\j,,dt < K{T), y T > 0, (1.16) 

te[o,T] Jo 

lim [ (\p- p\^ + p\u\^) {x,t)dx = 0, (1.17) 

t— >+oo J ^ ' 

lim / |Vu|2(a;,i)dT = 0, (1.18) 
where K{T) is a positive constant dependent on p_^, pi, \\{po — p, uo)||//3 and T. 

Remark 1.1. For example, we can choose that P = Ap'^ and A(p) = cp^ with 7 > 1 and /3 > 1, where A 
and c are two positive constants. Also, we can choose that A is a non-negative constant. 

Remark 1.2. Considering the case that the space domain C is bounded and A = p'', /? > 3, Vaigant- 
Kazhikhov ([H], Theorems 1-2) obtained the global existence of strong solutions when the initial data 
are large and the initial density is boimded from zero. In this paper, since the initial energy is small, we 
can use the similar argument as that in the case A =constant O [6l [7l [8] to obtain some good a priori 
estimates of the solution, and obtained the global existence of classical solutions when the space domain 
is R^, iV = 2 or 3, the initial density may vanish in an open set and /3 > 1. 

The proof of Theorem 11.11 bases on the derivation of a priori estimates for the local solution. Specifi- 
cally, in Section[2l we fix a smooth, local in time solution for which < p < p and A1+A2 < 2{Co + Cf)^ , 
then obtain the estimate Ai + A2 < (Co + C/)^, and prove that the density remains in a compact subset 
of (0,p). Using the classical continuation method, we can close these estimates. 

Using the initial condition uq G H^, we can obtain pointwise bounds for F in Proposition 12.61 which 
is the key point of the a priori estimates. Because that the mass equation can be transformed to the 
following form, 

^^A{p{x{t),t)) + P{p{x{t), t)) - P{p) = -F{x{t),t), (1.19) 

where A satisfies that A(p) = and A'(p) = Ii!±MPl^ a curve x{t) satisfies x{t) = u{x{t), t), thus pointwise 
bounds for the density will therefore follow from pointwise bounds for F. 

In theorem 1 1.21 the constant e is independent of p^. Thus, we can obtain the global existence of weak 
solutions to (|l.ip - (|1.3p with the nonnegative initial density po > (the two-dimensional result can be 
found in [20j). 

Definition 1.1. We say that (p, u) is a weak solution of (|l.ip - (|1.3p . if p and u are suitably integrable 
and satisfy that 



t2 rt 



p(j)dx 

tl Jti 

for all times ^2 > ii > and aU e C^^iR^ x [ii,t2]). 



+ pu-V(t))dxdt (1.20) 



pwipdx 
t 



*2 rt2 

{pu ■ ipt + p(w • V"!/;) • u + PdiY^pydxdt 



tl Jti 

{pdkU^ dkip-' + {p + A)divudivV' — pfip}dxdt (1-21) 



for all times ^2 > > and aU ip e (Co^(K^ x [ti,t2]))^ 
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Concerning the pressure P, viscosity coefficients /i and A, we fix < p < p and assume that 

'^^O' Mm), a^^s! f-aiipe[o,p], 

P(0) = 0, P'{p) > 0, (1-22) 



ip - p) [P{p) ~ P{p)] > 0, p e [0, p) u (p, p] , 
Pi-) 

2m+A(-) 



PeC^([0, p]) or 9^^w.-) is a monotone function on [0, p]. 



Theorem 1.2. Assume that conditions {i.6')] -i f77S)) and \1.22]) hold. Then, for a given positive number 

M (not necessarily small) and pi € {p,p), there are positive numbers e and 9, such that, the Cauchy 
problem (pUP^dUSP with the initial data (po,Uo) cind external force f satisfying 

< inf po < suppo < pi, 

Co + Cf<e, (1.23) 
Mq < M, 

has a global weak solution {p,u) in the sense of lll.20\) - lil.21\) satisfying 

C^Mnfpo < p< p, a.e. (1.24) 
p - p, pw e C([0,oo);i/~i), 



< " >m-7lToo) +sup(||VF(.,i)|U. + ||Vu;(.,t)||i.) < C(r)(Co +Q)^ (1.26) 
sup / |Vupdx+ / / a\Vu\^dxds <C, (1.27) 



T 



\\F{;t)\\L^dt<C{T), (1.28) 



where a e (0, 1) when N ^2, ae {0, i] w;/ien iV = 3, r > and a' G (0, ^^^], 

+ J J (|Vwp + cr|(pu)t+div(pu(K)u)p+CT^|V'ap)da:dt 
< (Co + C/)^ (1.30) 



Um (\p- p\'^ + pIu]"^) {x,t)dx ^0. (1.31) 
In addition, in the case that inf pQ > 0, the term J a\u\^dxdt may he included on the left hand side of 

Remark 1.3. Considering the case that A =constant, HofF-Santos and HofF O [7l [8] obtained the 
existence of global weak solutions. In this paper, since the viscosity coefficient A is a function of the 
density p, we need a higher regularity condition Vuo G L^, and use some new methods to obtain a priori 
estimates of the solution. Using the initial condition Vmq G L"^, we can obtain pointwise bounds for F 
in Proposition 12.61 which is the key point of the a priori estimates. Using the compensated compactness 
method [Hj and the estimate ||F(-, i)||L=odt < C(T), we can obtain the strong limit of approximate 
densities {p''}, see Sectional 
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Then, we study the propagation of singularities in solutions obtained in Theorem 11.21 Under the 
regularity estimates of the solution in Theorem II. 2( and similar arguments as that in ^ I20j , we can 
obtain Theorems 1 1 . 3f[T77l and omit the details. 

In Theorem ll.31 we obtain that each point of determines a unique integral curve of the velocity field 
at the initial time t — 0, and that this system of integral curves defines a locally bi-Holder homeomorphism 
of any open subset onto its image fJ* at each time t > 0. From this Lagrangean structure, we can 
obtain that if there is a vacuum domain at the initial time, then the vacuum domain will exist for all 
time, and vanishes as time goes to infinity, see Theorem 11.51 Also, in Theorem 11.61 we obtain that, if 
the initial density has a limit at a point from a given side of a continuous hypersurface, then at each 
later time both the density and the divergence of the velocity have limits at the transported point from 
the corresponding side of the transported hypersurface, which is also a continuous manifold. If the limits 
from both sides exist, then the Rankine-Hugoniot conditions hold in a strict pointwise sense, showing 
that the jump in the (A + 2/i)divM is proportional to the jump in the pressure (Theorem 1 1.7|) . This leads 
to a derivation of an explicit representation for the strength of the jump in A(p) in non-vacuum domain. 



Theorem 1.3. Assume that the conditions of Theorem \1.2\ hold. 

(1) For each to > Q and xq E M.^ , there is a unique curve X{-,xo,to) e Ci((0,c5o);]R^)nC^T^([0,oo);R^), 

satisfying 

X{t;xo,to) = Xq + / u{X{s;xo,to),s)ds. (1-32) 

J to 

(2) Denote X(t,XQ) = X{t;xo,0). For each t > and any open set ft C M^, f7* = X{t, ■)n is open and 

the map Xq i — > X(t,xo) is a homeomorphism of onto f2*. 

(3) For any < ti,t2 < T, the map X{ti,y) — > X{t2,y) is Holder continuous from M.^ onto M.^ . 

Specifically, for any yi,y2 G M.^ , 

\X{t2, y2) ~ X{t2,yi)\ < exp(l - e-^(i+^))|A(ii, ys) - X{t,,y,)r^'''^"' . (1.33) 

(4) Let M CR'^ be a {N - l)-mamfold, where (3 e [0, 1). Then, for any t > 0, = X{t, ■)M is a 

Cf^' [N - l)-manifold, where /?' = (}e-^^^+^\ 

Theorem 1.4. Assume that the conditions of Theorem ] hold. Let V be a nonempty open set in M.^ . 
//essinfpo|y > P > 0, then there is a positive number p~ such that, 

p{-,t)\v' > (f, 

for all t > 0, where V* = X{t, ■)¥. 

Theorem 1.5. Assume that the conditions of Theorem \1.2\ hold. Let U be a nonempty open set in M^. 
Assume that pq\ij — 0. Then, 

p{-,t)\u^ =0, 

for all t > 0, where — X{t, ■)U . Furthermore, we have 

lim \{x e R^\p{x, t) = 0}| = 0. (1.34) 

t — ^oo 

Recall that the oscillation of g at x with respect to E is defined by (as in [5]) 

osc{g;x,E) = lim (esssup5|£;nB„(x) ~ ess\nig\EnBR(x)) , 

R — '0 

where x € E and g maps an open set E C into M. We shall say that g is continuous at an interior 
point X of E, if osc(g; x, E) = 0. 

Theorem 1.6. Assume that the conditions of Theorem hold. Let E C be open and xo G E. If 
osc{po; Xq, E) = 0, then osc{p{-,t); X{t,XQ), X{t, ■)E) = 0. In particular, if xq G E and po is continuous 
at xq, then p{-,t) is continuous at X{t,XQ). 



5 



Now, let be a C° (TV — l)-manifold in and Xq G A4. Then there is a neighborhood G of a;o 
which is the disjoint union G = {G H A4)U E+ U where E± are open and xq is a limit point of each. 
If osc(g; xo, E+) — 0, then the common value g{xa+^ t) is the one-sided limit of g at xq from the plus-side 
of A^, and similar for the one-sided limit g{xo—,t) from the minus-side oi Ai. If both of these limits 
exist, then the difference [ff(a;o)] := g{xQ+) — g{xo~') is the jump in g at xq with respect to ^A (see [S]). 
Then, we can obtain the following result about the propagation of singularities in solutions. 

Theorem 1.7. Let {p, u) as in Theorem ] 1.2[ A4 he a {N — l)-manifold and xq G Ai. 

(a) If po has a one-sided limit at xq from the plus-side ofM, then for eacht > 0, p{-,t) and divu{-,t) 
have one-sided limits at X{t,Xo) from the plus-side of the C° (N — \)-manifold X{t, ■)Ai corresponding 
to the choice E*^ = X{t,-)E+. The map 1 1-^ p{X{t,xo)-i-,t) is m ([0, oo)) n C^((0, oo)) and the map 
1 1-^ divu(X(t, a;o)+, i) is locally Holder continuous on (0, cx)). 

(b) If both one-sided limits po{xoi) of po at xo with respect to M exist, then for each t > 0, the jumps 
in P{p{-,t)) anddivu(-,t) at X{t, xq) satisfy the Rankine-Hugonoit condition 

[i2^l + X{p{X{t, x„),t)))diYu{X{t, Xo), t)] = [Pip{X{t, Xo), t))]. (1.35) 

(c) Furthermore, if pq{xq±) > p > 0, then the jump in A{p) satisfies the representation 

[AipiX{t, xo), t))] = exp ^ a{T, xo)dT^ [Mpoi^o))] (1-36) 

where a(t xn) - IPM^jt^^o)^*))] 
wneie ayL,xo) - [A{p{x(t,xo),t))] ■ 

Remark 1.4. Using similar arguments as that in [20] . we also can show that the condition of p ^constant 
will induce a singularity of the system at vacuum in the following two aspects: 1) considering the special 
case where two fluid regions initially separated by a vacuum region, the solution we obtained is a non- 
physical weak solution in which separate kinetic energies of the two fluids need not to be conserved; 2) 
smooth solutions for the spherically symmetric system will blowup when the initial density is compactly 
supported. Thus, the viscosity coefficient p plays a key role in the Navier-Stokes equations. 

We now briefly review some previous works about the Navier-Stokes equations with density-dependent 
viscosity coefficients. For the free boundary problem of one-dimensional or spherically symmetric isen- 
tropic fluids, there are many works, please see [HI UHl HH [IZl HH US] and the references cited therein. 
Under a special condition between p and A, A = 2pp' — 2p, there are some existence results of global 
weak solutions for the system with the Korteweg stress tensor or the additional quadratic friction term, 
see [HE]. Also see Lionsfll] for multidimensional isentropic fluids. 

We should mention that the methods introduced by Hoff in 7J and Vaigant-Kazhikhov in [16 will 
play a crucial role in our proof here. 



2 Global existence 

Standard local existence results now apply to show that there is a smooth local solution (p, u) to (jl.ip - 
p.3|) . defined up to a positive time Tq, such that 



p, ue Ci(R^ X [0, To]) with p > for aU t e [0, To] 

and 

(p ~p,u) G C{[0, To],H') n C'{[0, To],H'). 

(See for example Matsumura-NishidajT4] and Nash 15|.) Let [0,T*) be the maximal existence interval of 
the above solution to (ll.l|) - (ll.3|) . 

In this section, we derive some a priori estimates for the local smooth solution of the system l|l.ip - (|1.3p . 

Claim 1: For any T > 0, if {p,u) satisfies 

0<p<p (2.1) 
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and 

Ai+A2<2{Co + Cf)^, yte [0,T]n [0,T*), (2.2) 

where 6 e (0, 1), 

Ai(T)= sup a / |Vupdx+ / / ap\u\'^dxdt 

te(o.T]n(().T*) J Jo J 



tG(0,T]n(0,T*) 

and 



p\u\^dx+ / / a^\Vu\^dxdt, T A T* = min{T, T*}, 



then we have 

ip^ < p < p Ai+A2< (Co + V i e [0, T] n [0, T*). (2.3) 

In this paper, we assume that e < 1. 

We can rewrite the momentum equation in the form, 

pu^ = djF + pdkw^^^ + pP. (2.4) 

Stated differently, the decomposition (|2.4p imphcs that 

AF = div(pu- p/). (2.5) 

Similarly, we have 

pAw;^'^- = afc(pu^) - a,(pit^-) + d,{pf) ~ dkipf)- (2.6) 

Thus estimates for pii, immediately imply hounds for WF and Vw. These three relations (|2.4p - (|2.6p 
will play the important role in this section. 

From now on, the constant C (or C{T)) will be independent of p^. 

Proposition 2.1. There is a positive constant C — C{p) independent of p^, such that if{p, u) is a smooth 
solution of U.l\) - U.3\) satisfying i2.1 \) -i2. ^) . then 



sup 

*e[o,T]n[o,T'; 



Ip\u\' + G{p) 





pT/\T' p 


dx + 


Jo J 



\Vu\^dxdt < C(Co + Cf). (2.7) 



Proof. Using the energy estimate, we can easily obtain (j2.7p . and omit the details. □ 

The following lemma contains preliminary versions of bounds for Vu and pit. 

Lemma 2.1. If{p,u) is a smooth solution of U.1\) ~ I[T7S\) as in Provosition \2.1[ then there is a constant 
C — C{p) independent of p^, such that 



t£(0,T]n(0,T') 



where Oi = Jq^"^ J crjVupdxdi, and 







1 \S/u\'^dx + 


Jo J 



ap\u\^dxdt <C {Co + Cf + 0i) , (2.. 



sup / plupdT + / / (7^ ( |Vu| 

ie(o,T]n(o,T') 



-— divit 
Dt 



dxdt 



< C{Co + Cf + Ai{T)) + C J J {\uf + \\'uf)dxdt. 



(2.9) 



Proof. In [20] (Lemma 2.1), we obtain this lemma in R^. Using the similar argument as that in 20J 
(Lemma 2.1) and [7J (Lemma 2.1), we can easily obtain this lemma in R'^ and omit the details. □ 

The following lemmas will be applied to bound the higher order terms occurring on the right hand 
sides of (HSl-lMl). 
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Lemma 2.2. If{p,u) is a smooth solution of lll.l\} - l[773\} as in Provosition \2. 1[ then there is a constant 
C — C{p) independent of p^, such that, 



2N-Np + 2p Np-2N 



[2,oo), N = 2, 



\\u\\l^<Cp\\u\\l^ IIV^II^,^' 'P^|[2i(J/' N = 3. ^'^■^^^ 

II IID ^ ^ 1^ ^ 2iV-Wp + 2p ,, 2W-Wp + 2p f [2,Oo), N = 2, 

(2.n) 

||Vu||lp <Cp(||F||iP + ||w||iP + ||P-P(p)||ip), (l,oo), (2.12) 
IIV^^IIlp + ||V?«||lp < Cp{\\pu\\L^ + II/IIlp), p e (1,(X3), (2.13) 

2Af- JVp + 2p 2iV-iVp + 2p 



\\F\\l. + \\w\\l. < Cp{\\pu\\^.'^ {\\Vu\\^, +\\P-P{p)h2 ) 

[2,cx3), / 
[2,6], N = 3 



WuU2 + \\fU2 + \\p-p{p)U2), pe{ g'^)' (2.14) 



Also, forO<ti<t2<T,p>2ands> 0, 

t2 pt2 



fa'\p~p\Pdx +C-^ f ' f a'\p- p\Pdxds 



(nl7{t2)Vti n pt2 n \ 

J J a'~^\p- p\Pdxds + J j a'\F\Pdxds\ , (2.15) 

j^^ j a'\p-p\Pdxds<c(^jJ j a'\FY'dxds + Co + C^ . (2.16) 

Proof. Using the similar argument as that in [20j (Lemma 2.2) and [7] (Lemma 2.3), we can easily obtain 
this lemma and omit the details. □ 

To bound the higher order term (t|VuP occurring on the right hand sides of (|2.8p in R^, we need to 
obtain the estimate of near t — Q. 

Lemma 2.3. If N — 3, uq G , {p,u) is a smooth solution of i f J. jp - f77g)) as in Proposition \2.1l then 
there is a positive constant Ti independent of p^ , such that 



sup / \Vu\^dx+ / p\u\^dxdt <C{1 + Mg). (2.17) 

l,TiAT]n[0,T*) JR3 Jo Jr^ 



sup / \Vu\'^dx+ / 

tG[0,TiAT]n[0,T*) JR3 J{ 

Proof. Using a similar argument as that in the proof of (|2.8p . we have 

/ |Vupdx+ / / p\u\^dxdt<C{Co + Cf+Mg) + c[ [ \Vu\^dxds. 

JR3 Jo JR3 Jo Jm.3 

From (|TT^ and ^JE\i . we have 

f f \Wu\^dxds<C + C f f {\F\^ + \w\^)dxds. 

Jo J«.^ Jo JwL^ 

From and ^JQ, we obtain 

/ {\F\^ + Iw]"") dx 

Jm3 

< c( [ iFl^dxY ( [ iVFpdxV+cf/ Iwl^dxY ( [ \\/w\^dx\ 
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< c(^Jj\Wu\' + \p-p\')dxy (^J^^{p\u\' + \f\')dxy . (2.18) 



Thus, from Proposition [2TT1 we have 



[ \Vu\^dx+ [ [ p\u\^dxdt 



< C(l + Mg) + Ct sup |iVu(-, s)\\l2, te [0, 1]. 
se[o,t] 

Thus, when Ti ~ minj gp^i^x^'^M yi 1 1}: we can easily obtain (|2.17p . □ 
Now, we apply the estimates of Lemma 12.21 to close the bounds in Lemma 12.11 



Proposition 2.2. // (p, u) is a smooth solution of lll.l\} ~ l[T7S\} as in Proposition \2.1l and s is small 
enough, then we have 



/^T/\T' p 
{a\S/u\'^ p\u\^)dx+ {ap\u\^ + a^\Vu\^)dxdt <{Co + Cff. (2.19) 

. . V - . "'0 J 

Proof. Since we obtain this proposition in in [20 (Proposition 2.2), then we only prove this proposition 
in M.^ in this paper. 

From Proposition 12.11 and Lemmas I2.1H2.21 we have 

LHS of (EUD < C(Co + C/) + C / / {a\Vu\^ + a^\u\^ + (T^\Vu\^) dxds. (2.20) 

From (|2.12p . we get 

/ / <T^\Vu\Uxds< / [\F\^ + \w\^ + \P - P{p)\^] dxds. (2.21) 

Jo Jo JB3 

From (1121), 113), (PTT|) . (jTTU)) and 11331) -(EUl) we obtain 

/ / {\F\'^ + \w\'^)dxds 

Jo Jr3 

/^^^ f / li^I'da;) ' ( [ \\/F\^dx) \( [ \w\^dx) ' ( [ \^w\^dx 

Jo \JR3 / VJR3 / VJR3 / \JR3 



j-TAT 

< c 



ds 



< C sup (/ ai\Vu\^ + \p- p\^)dx [ a^p\u\^ + \f\^)dx\' [ [ a {p\u\^ + \f\^) dxds 

te[o,T]n[o,T*)LjR3 Jr3 J Jo Jr^ 

< C{Co + Cf)^' + C{Co + Cff, (2.22) 

j-TAT* j- / j-TAT' j- \ 

/ / CT^Ip-pl^rfxds < C / / cr^li^l^dxds + Co + C/ 
Jo Jb3 yJo Jr3 j 

< C{Co + Cf)^' + C{Co + Cf), (2.23) 

j-TAT* j- j-TAT* 

/ / a^\u\'^dxds < C(Co + C/)3 / ct^ ^||Vw||i2 + HVwHis) ds 
Jo Jr3 Jo 

< CAUCo + Cf)i +CAi{Cn + Cf)i. (2.24) 

From (f23T|) - (12:241) ■ we have 

TAT* r 

cr3(|M|4 + iVM^fixds < C(Co + C/)2» + C(Co + C/). (2.25) 
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Similarly, we get 



TAT* 

3, 



TiATAT* 
TAT* 



a\Vu\ dxds 

pl AJ (■ 

< / / (a^lVul^ + \\'u\^)dxds 

Jt^ /\T/\T* Jr^ 

/ / {a^\Wu\-^ + \Wu\^) dxds 

JTif\T/\T' Jr3 



'TiATaT 

;.TaT 

< C(Ti) 



TiATAT* 



< CiA-QiCo + CfY" + C{M,){Co + Cf), (2.26) 



<7\\Iu\^dxds 



TiATAT* 
r-Ti aTaT 



< C{Co + Cf)+ a{ i\Vu\^ + \p~p\^)dx) / {p\u\^ + \f\^)dx) ds 

Jo VJr3 / VJr3 / 



< 



3 / ^TAT* ^ 

C(M,)(Co + Q) + C(Af,)(Co + C/)3A? + ( a\\^u\\l,ds 

TAT* i-TAT* 



pi AJ /• 

sup "'II'^'^IIl^ 

i,TiAT]n[0,T*) Jo "'0 



\te[o,' 

< C(M,)(Co + C/) + C(Mg)(Co + C/)^(i+^) + C(M,)(Co + C/)3+''. (2.27) 
Then, from ([121), and (p:^ - (p:?7)) . we obtain 

LHS of (ITTgil 

< C(M,)(Co + Q)i^2e^f(i+e)A(i+e) 

< (Co + C/)^ (2.28) 

when 

C(M,)£(i-«)^'''^(3-ie)Ai < (2.29) 

□ 

Then, we consider the Holder continuity of u in the following lemma. 
Lemma 2.4. Let a e (0, 1) when TV = 2, a e (0, i] iw/ien iV = 3. When t e (0, T] n (0, T*), we have 

/ iV-2 + 2e, 4-]V-2a 4-JV-2c, 1-q\ 

<u{-,t)>'^<C[\\pu\\^,^ (l|Vu||^, ^ +{Co + Cf)^^) + \\Vu\\L2 + {Co + Cf) — j. (2.30) 

Proof. Let p = From (|2.12p . (|2.14p and Sobolev's embedding theorem, we have 

< u{-,t) >" 

< C||Vu||lp 

< C(|1F||lp + ||w||lp + ||P-P(p)||lp) 

/ Np-2N 2N -Np + 2p 2N -Np + 2p 

< C i^WpuW^.'^ (||Vw||^. +\\P-Pip)\\L^ ) + ||V^i|U2 + ||/|U2 + ||P-P(p)|U2 

+c\\p^p\\Up-p0 

/ Wp-2iV 2JV-iVp+2p 2JV-]Vp + 2p l\ 

< cI^IIHIl^'" (llV^lli. +iCo + Cf)^^) + \\VuU2 + iCo + Cf)-^y 

□ 
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Proposition 2.3. If uq € , {p,u) is a smooth solution of f J. j)j -i [77g)) as in Provosition \2. 1\. then we 
have 

sup \Wu\^dx+ p\ufdxdt<C{Mq). (2.31) 

ie[o,T]n[o,T*) J Jo J 

Proof. Using a similar argument as that in the proof of (j2.8[) . we have 

sup / iVwpdx + / / p\u\'^dxdt < C(Co + Cf + Mg)+C / \Vu\^dxds. 

te[o,T]n[o,T*) J Jo J Jo J 

From Lemma [^751 and (|2.26p . we can easily obtain (|2.3ip . □ 

Proposition 2.4. If uq E H^, {p,u) is a smooth solution of f J. jp ~i fI3)) as in Provosition \2.1\. then we 
have 

sup a p\u\^dx+ I (j\S/u\^dxdt <C{Mq). (2.32) 

tG(0,T]n(0,T*) J Jo J 

Proof. Since we obtain this proposition in in [20 (Proposition 2.4), then we only prove this proposition 
in R'^ in this paper. 

Using a similar argument as that in the proof of (|2.9p . from (|2.3ip . we have 

/. /.TAT* j- j-TAT' j- 

sup cr / p\u\^dx+ / cr\\/u\^dxdt<C{Mg) + C / ct + IVul"*) dxds. 

te(0,T]n(0,T*) JR3 Jo Jo Jr^ 

Without loss of generality, assume that T > 1. From (|2.25p . we get 

sup aj p\u\'^dx+ / a\Vu\'^dxdt<C{Mg) + C / cr + |Vm|'') da;ds. 

tG(o,T]n(o,T*) JR3 Jo Jr3 Jo Jr^ 

From (1121), (piTll - fTT^ and (P3T1) . we have 

sup a p\u\'^dx+ / a\Vu\'^dxdt<C{Mg) + C / crdFl"* + |u;|'')dxds. 

tG(0,T]n(0,T') JR3 Jo JR3 Jo JR3 

From dlJl), (1TT3| . (l2l0l) and (I23T1) . we obtain 

r-lAT* 



Jo JR3 



< / ^ f / l-Fpda;) ' ( [ iVFpdx) ' ds + C [ a( f \w\^dx] ' f / iVwpda;^ ' ds 

Jo \JR3 / \JR3 / Jo 

r-lAT 



< ^ a(^Jj\Vu\^ + \p-p\^)dxy (^J^^{p\u\' + \f\^)dxy ds 



< C{Mg)+C{Mq) sup a--\\^pu\\L^. (2.33) 

tG(0,T]n(0,T*) 

Using Young's inequality, we can finish the proof of this proposition. □ 
Lemma 2.5. For any p G [2, oo) when N — 2, p E [2, 6] when N ~ 3, we have 

2JV- Afp + 2p Np-2N 

MLp<Cp\\y^u\\^, \\\7u\\^,'- +Cp\\\/u\\l2. (2.34) 



Proof. Since 



p J \ufdx< J p\ufdx+ (^j Ip-ppd^y (^j \u\'^di^ 



applying (|2.10p . we get 

\\u\\l^<C\\^pu\\l.+C\\Vu\\l.. 
From (|2.10p . we can immediately obtain (|2.34p . □ 
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Lemma 2.6. For any q € (0, 2) when iV = 2, 5 S (1, |) when N = 3, we have 



TAT* 

' aP'-^p\u\^+'^dxds < C{Mg). (2.35) 







Proof. Since we obtain this lemma in in [20j (lemma 2.5), then we only prove this lemma in MP in 
this paper. 

Using Holder's inequality, ((OT|) . ((O^ and ((Oi)) with p = 6, we have 



pT/\T P 

/ / a'-^p\u\^+'^dxds 

Jq JR3 

< C / fj 4 WVp^WJ \\u\\L6ds 











L2 II V u,||^2 

3q 4-3g 
^TAT* \ * / /-TAT* 



< C(Afg). 



□ 



Proposition 2.5. // uq G -ff^, (Pi w) is a smooth solution of m.l\) - nT3\) as in Proposition\2.4\ and 



< wt^^ , Vp£ [0,p], (2.36) 







^ pTat' !■ 

sup <jP' p\u\^+''dx+ aP'\u\'^\Vu\'^dxdt <C{Mg). (2.37) 

t6(0,T]n(0,T') J Jo J 

Proof. In [5D] (Proposition 2.5), we obtain this proposition in R^. Using the similar argument as that in 
[20] (Proposition 2.5), we can easily obtain this proposition in and omit the details. □ 

Proposition 2.6. // / G Lf'L/^'^, {p,u) is a smooth solution of U.l]) - ^^) as in Provosition \2.5[ then 
we have 

2(2 + g-W) 2iV + JVq 

\\f\\l^ + \\w\\l^ < c{\\Wu\\l2 + Hp- p|U2)^ww(||pu|U2+, + 11/11^2+,)™^ (2.38) 

and 

tat* 

(llFllioc, + |lu,||ioo)ds< C(M,)(Co + C/)'*+^.V«J(l+T). (2.39) 
Proof From (P??)) . ((TTO)) . ([^1^ . (|07)) and the Gahardo-Nirenberg inequahty, we have 

\\F\\l^ 

2(2 + g-iV) 2JV + Afg 

< c\\F\\ir^+''^ wvFwitit"' 

2(2 + g- Af) 2iV + iVg 

< C(||VM||i2 + \\p~p\\l2) 4 + 2, + iV, (||pw||i2 + , + ||/||i2 + ,)*+ = ' + ~<'. 

<.tat* 

/ \\F\\L^ds 
Jo 

r'^^'^ ^1 g 2(2 + g-JV) pi 2JV + Wi; 

< C{Mq) / (cr 2 (Co + C/) 2 ) 4+25+iv, (cr 2+, ) 4+2,+N, (Is 

Jo 

< < C(M,)(Co + C/)^ww(l + r). 

Similarly, we can obtain the same estimates for w. □ 



and 
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Then, we derive a priori pointwise bounds for the density p. 

Proposition 2.7. Given numbers < < p < pi < p2 < p, there is an e > such that, if {p,u) is a 
smooth solution of 11.1\} - ![T73\) with Co + C/ < e o,n-d p^ < po < pi, then 

\p^<P<P2. (x, i) e X {[0, T] n [0, T*)}, (2.40) 
for any T > 0. Furthermore, Claim 1 and the estimates in Propositions 1 2. I\i2. 6\ hold for any T > 0. 



Proof At first, we prove that if ((2?T1) and ((2?2)) hold, then estimate (|2^ holds. 
We fix a curve x(t) satisfying x — u{x{t), t) and a;(0) = x. From (|1.19p . we have 

^A(p(x(t), t)) + P{p{x{t), t)) - P{p) = -F{x{t),t), (2.41) 

where A satisfies that A(p) = and A'{p) = ^^ii^. 

(I) For the small time, we estimate the pointwise bounds of the density as follows. From (|2.ip and 
((O^ . we have, for all t e [0, 1], 

\A{p{x{t), t)) - A{po{x))\ < C{Mg){Co + Q)^S^S^ + Ct. 

When 

2C{Mg)e^m^ < A(pi + -(p2 - Pi)) - A(pi), (2.42) 

we get 

A{p{x{t),t)) < A(pi + i(p2 - pi)), t e [0, r], 

and 

P<Pi + liP2-Pi), (x,t) gM^ X [0,t], (2.43) 
where r — min{l, ^[A{pi + ^{p2 — Pi)) — A{pi)]}. Similarly, since 

Hp,)-Hj:P,)> I -^ds = 2fAn-, 



then, if 



we get 



2C(M,)e^SOTr < 2/x In | < A{p^ ) - A(|p J, (2.44) 

P>Ipv (a;,t)eR^x [0,ri], (2.45) 

where n = minjr, ^ In |}. 

(II) For the large time t > ri, we estimate the pointwise bounds of density as follows. From (|2.7p . 
(|2l^ . ((OTI) . (|238| and (EHI), we have 

^^^MMM)+p(,(,(,),,))_p(^).0.(i), (2.46) 

where 

|05(<)| <C(Ti,M,)(Co + C/)WW, f>ri. 
Now, we apply a standard maximum principle argument to estimate the upper bounds of density. Let 

to ^ snp{t e (r,T] n {t,T*)\A{p{x{s), s)) < A{p2), for aU s £ [0,t]}. 

If to < I" and to < T* , we have 

A(p(a;(to),to)) = A(p2), 
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dAipix{t),t)) 



dt 

and 



>0, 



p{x{to),to) = p2. 

From (|2.46p . we have 

O^ito) > Pip2) ~ P{p). 

On the other hand, when 

C(ti,M,)£^+55+w < P(p2) - P(p), (2.47) 

we have 

O^ito) < P{P2) - P{P). 

It is a contradiction. Thus, we have 

p < p2, (a;, e X {[0, T] n [0, T*)}. (2.48) 

Similarly, let 

<i = sup{i e (t,T] n (T,r*)|A(p(x(s),s)) > A(ip^), for aU s e [0,i]}. 
If ti < T and ii <T*, we have 

A(p(a;(^i),ii))=A(ip^), 
dA(p(a;(t),i)) 



dt 

and 



<0, 
t=ti 



p(a;(ii),ii) = ^p^. 
From (|2.46p . we have 

O^ih) < P{l:p^) - Pip) < max P(s)-P(p). 

On the other hand, when 

a(2 + g-JV) 

C(ti, M,)£ < P{p) - max P(s), 

sg[o4p] 



05(^1) > max Pis) -Pip). 
se[o,ip] 



we have 

It is a contradiction. Thus, we have 

P>Ipv (a;,i)eR^x{[0,T]n[0,T*)}. 

Using the classical continuation method, (|2.19[) and (|2.40p . we can finish the proof of this proposition. □ 

From now on, the constant K iKiT)) will depend on p^ (and T). 
Lemma 2.7. For any T > 0, we have 

T/\T' 

||Vp(-,i)||i2+, + / 

te[o,T]n[o/ 



sup ||Vp(-,i)||i2+, + / (||Vu||l- + ||Au||i2+,)rfi < if(T), 
,Tln[o,T*) Jo 



(2.49) 



sup ||Vp(-,<)|U2 <X(r). (2.50) 

tG[0,T]n[0,T') 
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Proof. From (|1.19p . we have 

dtA{p{x, t))+u-WA + P{p(x, t)) - P{p) = -F{x, t), (2.51) 
where A satisfies that A(p) = and A'(p) — 1tl^iEl_ By the simple computation, we have 

||VA(t)f^t\ < \\vm\\ltl.+K f (|1VA|1^+^, + |lVF|U.+,|lVA|li+^, + |1 V^IU^ |1 VA|l^t\) ds. 



Using the Fourier analysis methods, one can obtain the following estimate. 

||Vu||loo < C|1m||l2 +C(||Vu||bo^^ + l)log(e+ ||Au||i2+,). (2.52) 
(For the convenience of reader's reading, we also give the proof in Appendix (|5.ip .) From (|1.12p . we have 

||Aw||i2+, < i^(||V^^|U2+, + ||F|U^||VA|U2+, + IIVAIU2+, + \\Vw\\l2+.), 

and 

< C{\\F\\l^ + \\p - PWl^ + MM- 

Thus, we have 

||VA(0||i+«, < ||VA(0)||^+^,+i^ f \\WF\\L2+,\\VK\\^+l,ds + K f A\og{\\V + e)\\V kf+^ds, 

Jo Jo 

and 

sup ||VA(t)|li2+, < (if+|lVA(0)|li2+, +if / |lVf|li2+,ds ) 

iG[0,T]n[0,T*) Y Jo J 

where yl= (||i^||L- + + 1) log (||i^||L- + ||u)|1loo + \\VF\\l2+, + \\Vw\\l2+, + e). From ([13), 

((2?3T|l . (|2?35| . ((07l) . ((OSl) and ((2^ . we can immediately obtain (|2^ . Similarly, we can obtain 
(1230)1 . □ 

Lemma 2.8. If po — p E and uq G , then for any T > 0, we have 



sup / p\u\'^dx + [ I f |VmP + l^divup^ dxdt < K, 
rin[o,T*)J Jo J \ Dt J 



te[o,T]n[o,T 



(2.53) 



sup {\\u\\l^. + \\u\\h2 + \\VF\\l2 + \\Vw\\l2) <K{T). (2.54) 

tG[0,T]n[0,T*) 

Proof. Using the similar argument as that in the proof of Proposition 12.41 we can obtain (|2.53p . From 
(mni), we get 

sup (||Vi^||i2 + ||Vii;||L2) < if. (2.55) 
te[o,T]n[o,T*) 

From jlJl), (frT2l) . (12311) . (IM) . (IM) . (l230l) . (l235l) . we have 

\\u{;t)\\H'^- < Ki\\u\\L2 + \\VF\\L2 + \\FVp\\L2 + \\Vp\\L2 + \\Vw\\L2) 

2 + q-N N 

< K{\\u\\l2 + ||VF||i2 + llV^^ll^r ||Vp|U2+, + ||Vp|U2 + \\\/w\\l2) 

< K{T), t e [o,r] n [o,T*). 

Then, using Sobolev's embedding theorem, we can finish this proof. □ 
Lemma 2.9. For any T > 0, we have 

sup \\pi;t)-p\\H2 <KiT). (2.56) 
te[o,T]n[o,T*) 
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Proof. From (|2.5ip and the simple computation, we have 

\\Ait)fH2<\\A{0)fH2+K f {\\F\\H2\\A\\m + {l + \\Vu\\L^ + \\u\\H2)\\AfH.)ds. (2.57) 



From (P3)l . ([QT]) . ([2301), (P^O)) . (P35)) . the Gahardo-Nirenberg inequahty and Sobolev's embedding 
theorem, we get 

\\F\\h- < K{\\F\\L2 + \\Vpu\\L^ + \\pVu\\L2 + \\Vf\\L2 + \\f -Vph^) 

< i^(r)(i + ||Vp||L3|iu|ii6 + ||Vw|U2) 

/ 6- W JV 3-iV JV \ 

< i^(r)(i + ||vp||^i \\v'p\\i4u\\J liv^ill^^ + liv^iiu^j 

< K(T)il + \\A\\H2{l + \\Vu\\L2) + \\Vu\\L2). (2.58) 
Thus, from (f2^ . ([2?53l) . (l237|) - (|2?58l) and Gronwall's inequahty, we have 

\\A{t)\\l,<K{T)+K{T) f {l + \\Wu\\L^ + \\Wu\\L2)\\A\\l2ds 

Jo 

and 

WAml^KKiT). 

Using (|2.40|) and (|2.50[) . we can immediately obtain (|2.56|) . □ 
Lemma 2.10. For any T > 0, iwe /lawe 

\\u\\jj,dt<K{T). (2.59) 



Proo/. From jM]), (l233l) . ([236]) and ([2381) . we have 



(||F||j^2 + ||u;||H2)'dt<if(T). (2.60) 
From jMl), (frT2ll . (|236l) and ((2J0| . we have 

/ |k|Pff3dt<i^ / (||M|U2 + (l + l|F||H2)(l + ||p-p||H2) + ||w||H2)2di <X(r). 

./o Jo 

□ 

Proposition 2.8. For any T > 0, we have 

rTAT" 



te[o,T] 



|Vupdx+ / \V^u\'^dxdt < K{T), (2.61) 



TAT* 

2 



sup {\\ip-p,u)\\„s + \\{pt,ut)\\H2)+ \\u\\j,,ds<KiT). (2.62) 

te[o,T]n[o,T*) Jo 

Proof. Taking the operator Vc)t + Vdiv(u-) in (jl.ip 9. multiplying by Vfi and integrating, we obtain 
p\\/u\'^dx 

= polVuopdx + J J {~\7pdtu'^u — \/{puj)dju\7u 

-Vii^VldjPt + dividjPu)] - pAu^[Aui + div(uAu^')] 

-Au^[djdt{{X + /i)divu) + div(udj{{X + ^)divu))] - Au^[{pf^)t + div{upf^)]}dxds 
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1=1 

Since po — p ^ H^, pa G [p^, pi] and uq G H^, we get 



(2.63) 



(2.64) 



Using (|2.40p . (|2.53P " (|2.54p . (|2.56p . (|2.59p . the integration by parts and Holder's inequality, we have 
J2 — — / pdfu'^udxds 



/^A-it + V((/i + A)divM) - VP - p/ 
V P 



Vudxds 



< K{T) J J \Vp\\Vu\{\V^u\\Vu\ + |Vu||Vpr + iVu^lVpl + \Vp\' + |Vp||V^u| + 

+ |V2p| + |Vu||Vp| + iV^wl + |V2p||Vu| + |Vp||-^divu| + |Au| + |V^divu| + \ft\)dxds 



Dt 



D 



< K{T) |||Vtt||2.||V^u|||,(||V^w|U2 + ||V— div«|U2 + ||/,|U2)||p-p||H2 

D 



3-N 

MVu\\J WV^ii 



L2 



+ 1) (IIVpll?,, + 1) + \\Vp\\h4^\\hs + 



Dt 



-divM 



IIVpll 



L2 



< K{T) + ^J^ (\\^'n\\h + llV^div^^lli, ) ds, 



h = - 



ds 
(2.65) 



j J V {pUj)djvSIudxds 

< c[ ||Vu||i4 (||Vp||l2||w||l- + ||Vu||i2)da;ds 
Jo 



(2.66) 



J4 = J Au^[djPt+diY{djPu)]dxds 

= - J [djAii^P'pt + dkAu^djPu'']dxds 

" i) / [^''^^^^"^J^^^ ~ i9A^'u'')P + Pdj {dk^ii^ u^)]dxds 

< J^l^J^ J (iV^ul + |Vp||Vu|)^da;dsy (^J^ J \W^u\^dxdsy 

J J pAu^[Aul + div{uAu^)]dxds 
J n[d^Au^d^u{ + Au'u ■ VAu^]dxds 

J ^hlV^Mp - d,Au^u''dkd,u^ - d^Aii^d^u'^dkU^ + Au^u ■ VAu^jdxds 



(2.67) 
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Jq — ~ i I Au^[djdt{{X + n)divu) + div{udj{{X + fi)divu))]dxds 



{djAu^[dt{{X + fJ.)divu) + div(ii(A + ^)divu)] 

+Au-'div(9ju(A + fi)diYu)}dxds 
t r 

djAu'[dt{{\ + ^)d\vu) + vf'X' dkpdiYu + (A + ii)u'^ dkdivu]dxds + Oq 
D 







djAu-'[{X + /i)— divu + X'ptdiYu + u~ X' dkpdivu]dxds + Oq 



I I djA{dtu' + u ■ Vu^){X + p)^dYvudxds + Oq 
Jo J Dt 

- [ + f^)\^§-^dWu\'dxds + K{T) + ^J^ W^Mhds, (2.69) 

J7 = -J^ J Aii^ [{pf)t + div{upf)]dxds 
J Aii^ipfi + pu-S/f]dxds 



- -7: [ [ P\^^u\'^dxds + K{T), (2.70) 
10 Jo J 

where Og denotes any term dominated by C /g /(|V/9| |Vup + |Vw||V^u|)|V^u|dxds and O7 denotes 

any term dominated by C J* J {\^divu\\V p\ + |Vu||V2-u|)|V-^divM| + \Vu\\V^u\\V p\\^divu\dxds, t e 
[0, T] n [0, T*). From ^^-^JO^, we immediately obtain Using similar arguments as that in the 

proof of Lemmas l2.9H2.10l we can easily get (|2.62p . □ 

Using the standard arguments based on the local existence results together with the estimates (|2.40p 
and (|2.62p . we can obtain that T* = 00. Since the uniqueness of the solution {p— p, u) S C([0, 00); H^) n 
C^([0, 00); i?^) is classical, we omit the detail. Thus, we finish the proof of the existence and uniqueness 
parts of Theorem 11.11 



3 Large time Behavior 

From (|2T6| . ((2T9)) and ([2^, we have 

J{\p- p\^ + \F\^)dxdt<C. 

From (|2.15p . we have 

r rN+2 

\p- p\^ix,t)dx < J \p- p\^ix,s)dx + C J I \F\^dxdT, 



(3.1) 
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where t € [iV + 1, iV + 2] and s G [N, N + 1], N > 1. Integrating it with s in [N, A'' + 1], we obtain 
sup / \p ~ pf{x, t)dx <C i\p-p\'^ + \F\^)dxdT. 

te[N+l,N+2]J Jn J 

Letting N oo, using (13. ip . we can easily obtain 



From (|2.25p . we can obtain 



From (|2.7p . we have 



hm / \p — p\ {x,t)dx — 0. 

t—>+oo 



J {lul"^ + \\7u\'^)dxds <C. 



J \Vu\'^dxds < C. 

Thus, for all e G (0, 1), there is a positive constant T^, such that for all r > Tg, we have 

{\u\'^ + \Wu\^ + |Vup + l/l^) dxds < e. 

For alH > Te + 2 and re [t-l,t-2], from (fTT]) . (fO]) and Holder's inequality, we get 
^1 



(3.2) 
(3.3) 

(3.4) 



< 



{x, t)dx + / |Mp [Ai|V-up + (A + ^)(divu)^] dxds 



1 



"1^") - 2^'^'"'^'^ ^ + M)diviiii • V|w|^ + pf ■ u\u\^]dxdt 



{x,T)dx + J J [Pdiv(| 
{x,T)dx + C (^J^ J iVupdxdsj (^J^ J \u\'^dxds^ 



-C 



J \Vufdxds^ ' (^J^ J \u\^dxds^ ' + (/ / l/l^^2;ds^ ' J \u\'^dxds^ 



< 



(a;, T)dx + Ce. 



Integrating it with r in [< — 1, t — 2], we obtain 

rt-i 



(x, t)dx < 



t-2 



(x, T)dxdT + Ce < Ce. 



Thus, we immediately obtain (jl.l7p . 
From ()2.34p . we have 

\\u\\l^ < C\\y^u\\L2+C\\Vu\\L^. 

From (|2.19p . (|3.2p - p.3p and p.Sp . we have that for all e G (0, 1), there is a positive constant T2c, such 
that for all r > T2e, we have 



(3.5) 



^ J {\uf + \Vu\^ + \Vu\^ + \u\^ + iViil^) dxds < e. (3.6) 
For all t > + 2 and t e — 1, t — 2], multiplying (jl.ip ^ by u, integrating it over x [r, i], we obtain 



p\u\^dxds 



{—ii ■ VP + p,Au ■ ii + V((A + /x)divu) ■ ii + pf ■ u) dxds 
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^J.. (3.7) 



Using the integration by parts and Holder's inequality, we have 

Ji - / [ ii- VPdxds 



J J divu(P - P{p))dxds 
j j \Vu\^dxds 



< C 

< Cei, (3.8) 



J2 = I I lJ,^u ■ iidxds 



I / \Vu\^{x,t)dx + !^ I |Vu|2(a;, T)dx ~ / / fidiU^ di{u''dkU^)dxds 



- ^2 J + I J \Vu\'^{x,T)dx + C J J \Vu\^dxds 

- "2 / + f / \VuWx,T)dx + Ce, (3.9) 

J3 = J y V((A + fi)divu) ■ iidxds 

= y" [(A + ^)|divMp] (x,t)da; + i y" [(A + /i)|divitp] (x,T)dx + J J ^X' pt\divu\'^ dxds 
t p 

(A + /i)divudiv(u • Vu)dxds 



< -\J [(A + Ai)|divMp] (x,t)da; + i J [(A + /i)|divMp] (x, T)dx + J \Vu\^dxds 

- + (2^'*)^2; + iy"[(A + /i)|divMp] (x,T)da; + Ce, (3.10) 

J4 = y y p/ • iidxds < C (^y y lup^xdsj < Cei (3.11) 
From jSJl), ((3Jl) - ([3lT|) . we obtain 

^ \Wu\^{x,t)dx + [ [ p\u\^dxds <Ce^^ +C [ |Vup(a;,r)dx 



2 

Integrating it with t in [< — 1, t — 2], we obtain 



^ y \Wu\^{x,t)dx <Ce^ +C J |VMp(x,T)dx < Cei 



Thus, we immediately obtain (|1.18p . 

Thus, we finish the proof of Theorem ll.il 
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4 Proof of Theorem 11.21 

Let jsix) be a standard mollifying kernel of width S. Define the approximate initial data (pq, Uq) by 

Pa = js *Pq + S, Uq = js * Uq. 

Assuming that similar smooth approximations have been constructed for functions P, f and A, we may 
then apply Theorem 11.11 to obtain a global smooth solution {p^,u^) of (|l.ip ~ (|1.3p with the initial data 
{Pq, Uo), satisfying the bound estimates of Propositions 12 . 1 j | 2 . 71 with constants independent of S. 
First, we obtain the strong limit of {u*}. From (|2.19p and (|2.30p . we have 



<u\-,t) >°'<C{t), t>T>0, 
where a G (0, 1) when N — 2, a G {0, i] when N = 3. From (|4.ip . we have 

1 



(4.1) 



u''ix,t) - 



Br{x) 



Br{x) 



u {y,t)dy 



< C(T)i?", t > T > 0. 



Taking R=l, from ((TTT|) and ^J^, we have 

l|w''||L~(R"x[r,oo)) < C{t). 



(4.2) 



Then, we need only to derive a modulus of Holder continuity in time. For all t2 > ti > r, from (|2.7p . 
(|TT9| . ((Oil) and g^l), we have 



y{x,t2)-u\x,ti)\ 

rt2 



< 



1 



\Br{x)\ J,, 
< CR-^\t2~ti 



Br(x) 



\ul{v,s)\dyd.s + C{T)R'' 

1 

\ul\^dyd^ +C{t)R°' 

1 

\u^\'^ + \u^ ■Vu^\'^dyds\ +C{t)R° 



< CR--\t2-ti\^ 

< C{T){R-^\t2-ti\i +R''). 
Choosing R = \t2 — ti \ "+2= , we have 

From the Ascoli-Arzela theorem, we have (extract a subsequence) 

-> w, uniformly on compact sets in X (0,oo) 



(4.3) 



(4.4) 



Second, we obtain the strong hmits of {F^} and {w^}. From (pi^ - ipH)) . and ((07)l . using 

similar arguments as that in the proof of (|4.ip - (|4.2p . we have 

< F\;t) +\\F'\\L^iR-^[r,T])+ < w'{;t) +\\w'\\L^^R.^[r^T]) < C(r, T), (4.5) 

where < t < t < T and a' e (0, ^1^*^^ ] . The simple computation implies that 



Ff = p\2^, + Xip'))lF' 



ds V2/i + A(s) 



^d_ / P{s)-P{p) 

<5« „,5 



divM 



-u' ■ VF' + (2/i + A(p''))divM' - (2/i + X{p'))^^u'^^JU.^ 



and 



(w^)^-^ = -u^ . W{w^)'''' + d.jui - dku^j - djufd^ui + dkufd.u 



(4.6) 
(4.7) 
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Then, from (ITT3| . ([2l9)) . (|2:25)l . gj) and gS]), we have 

II-P'/||l2(R«x[t,T]) + ||L2(jj,jVx[r,T]) < C(t,T), T > T > 0. 

Usmg a similar argument as that in the proof of (j4.3p . we obtain 
and (extract a subsequence) 



F" F, w ~^ w, uniformly on compact sets in M x (0,oo). 
Third, we obtain the strong limit of {p^}. From (|2.40p . we get (extract a subsequence) 



(4. 



(4.9) 



/^p, weak-* in L°°{R"). 

Let $(s) be an arbitrary continuous function on [0,p]. Then, we have that (extract a subsequence) ^{p^) 
converges weak-* in L°°(M.^). Denote the weak-* limit by 3>: 



From the definition of F, we have 
where 



$(p*)^$, weak-* in L°°(M^). 
divw = DF + Pq , 



1 



2p + A(p) 



, Poip) = ,^{p)iP{p)-P{p)). 



From (jl.ip . we have 
and 



dtplnp + div{p In pu) + Fpiy + pPo — 



9t(pln p) + div(plnpu) + Fpv + pPo = 0. 



Letting 'J = plnp — plnp> 0, we obtain 

+ +div(«'?i) + F(pF - pv) + Fp(t/ - P) + ^ - z?;^ = 0. 
with the initial condition 'I'|f=o — almost everywhere in M^. Let (/)(s) = sins. Since 



we get 



and 



0"(s) = ->i 5e[0,p], 

s p 



</.(p^) - m = ^'ip)ip' -p) + 20"(p + ap' pW - pf, e e [0, 1], 



lim||p*-p||i. <q|vi/|u,. 

o — !-U 



< C / |<?|vl/dx, 



Similarly, every function / G C^([0,p]) satisfies 

9{f-f{p))dx 

where g is any function such that the integrations exist. Then, when v G C^([0,p]), we have 

j F{pU~pi^)dx <C j \F\^dx 



and 



Fp{v — v)dx 



< C / \F\^dx. 



(4.10) 



(4.11) 



(4.12) 
(4.13) 

(4.14) 
(4.15) 
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When Po € C^i[0,p]), we have 



(pPo - pPo)dx 



< 



{pPo - pPo)dx 



p(Po - Po)dx 



< C / 



(4.16) 



When Pq is monotone function on [0,p\, using the Lemma 5 in [16J, we have 

'pPo>pPo. (4.17) 

From (|iTT|l - (|iT7)l . we obtain 

J "idx < J{1 + \F\)^dxds. 
Using (|2.39p and Gronwall's inequahty, we get 

* = 0, {t,x) e [0,T] X R^, 

and (extract a subsequence) 

p^ ~p^ p-p, strongly in L'=(R^ x [0, oo)), 

for all k e [2, oo). 

Thus, it is easy to show that the limit function (p, u) are indeed a weak solution of the system (jl.ip - 
(jl.SP . Using a similar argument as that in the proof of (|1.17p , we get (|1.31|) . This finishes the proof of 
Theorem [m □ 



5 Appendix 

It requires a dyadic decomposition of the Fourier space, so let us start by recalling the definition of the 
following operators of localization in Fourier space (see [3]): 

A,a^.F-i((^(2-«|C|)a), for 9 e Z, 

where Ta and a denote the Fourier transform of any function a. The function Lp is smooth, and satisfies 

supptp G |e e lei e [J J] I , E -^(^-^'i) - 1, V t G R\{0}. 
Let us note that if |j — j'| > 5, then s\rp^tf{2~H) n supp(/5(2^^ i) = 0. 

Definition 5.1. We denote by ^ the space of distributions, which is the completion of 5(R^), N > 2, 
by the following norm: 

dot ||r,sA;|| A II II 



all. 



Lemma 5.1 (p]). Denote B a hall of S.^ , and C a ring of M.^ . Assume that 1 < P2 < Pi < 00 cmd 
1 < 92 < 91 < OO- If the support of d is included in then 

II II ^ M\a\+N(-^ I 

// the support of d is included in 2^C, then 

< 2-'=*' sup ||a"a|lLPi. 

|q|=M 

Lemma 5.2. For any p > N and N > 2, there exists a positive constant Cp such that 

II V^IU^ < C\\uU. + Cpdl V^bo + 1) log (e + J . (5.1) 
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Proof. 

M-l 

\\Vu\\l^ < C\\u\\l2 + J2 llVAfcullioo = C\\u\\l^ + J2 l|VAfcw||ioo + ^ ||VAfeu||i= 

k>0 k>0 k>M 

< C\\u\\l^ + M\\Vu\\so + V 2*=(f-^)||AAfc«||iP 

oo ,00 ' * 



k>M 



< C\\u\\l2+M\\Vu\\so + j^\\Au\\lp. 

00 ,00 1 o "T: — 



1-2- 



Choosing M log + Hvlj l ^j^^'' +1 ^ ' ^® finish the proof. □ 
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